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Abstract. The inverse problem of diffraction theory in essence amounts to the reconstruc- 
tion of the atomic positions of a solid from its diffraction image. From a mathematical 
perspective, this is a notoriously difficult problem, even in the idealised situation of perfect 
diffraction from an infinite structure. 

Here, the problem is analysed via the autocorrelation measure of the underlying point 
set, where two point sets are called homometric when they share the same autocorrelation. 
For the class of mathematical quasicrystals within a given cut and project scheme, the 
homometry problem becomes equivalent to Matheron's covariogram problem, in the sense 
of determining the window from its covariogram. Although certain uniqueness results are 
known for convex windows, interesting examples of distinct homometric model sets already 
emerge in the plane. 

The uncertainty level increases in the presence of diffuse scattering. Already in one 
dimension, a mixed spectrum can be compatible with structures of different entropy. We 
expand on this example by constructing a family of mixed systems with fixed diffraction 
image but varying entropy. We also outline how this generalises to higher dimension. 



1. Introduction 

After 25 years of quasicrystal research, our understanding of the atomic structure of qua- 
sicrystalline alloys is still far from being complete [18]. The main reason for this is the 
difficult inverse problem of determining the structure at the atomic scale from the available 
information, which exists mainly in the form of diffraction intensities. Here, we discuss the 
non-uniqueness arising from homometric point sets, which is even present in the idealised 
situation of a perfect diffraction measurement from an infinite point set A C R d . 

First, we consider the situation where A is a mathematical quasicrystal or model set. A 
perfect diffraction image of A, as described by the positive diffraction measure^, uniquely 
determines its inverse Fourier transform, which is the autocorrelation (or Patterson) measure 
j A . The starting point is thus the (hypothetically complete) knowledge of j A , and the remain- 
ing task is then to determine A from this information. For a model set based on a known 
cut and project scheme, this amounts to determine the corresponding window W in internal 
space. 

Beyond pure point diffraction, we reconsider the known homometry between the binary 
Rudin-Shapiro sequence and the Bernoulli (or coin flipping) chain [10] . We introduce a new 
process, called 'Bernoullisation', which provides a continuous isospectral transition between 
these two extremal cases. This method generalises to arbitrary dimension and shows that 
even a perfect diffraction image (of mixed type) may not be able to distinguish structures of 
different entropy. 
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2. HOMOMETRY 



For finite point sets F C R , homometry is defined in terms of their difference sets F — F, 
taking into account multiplicities. Two finite point sets are called homometric when they 
share the same weighted difference set (which is a multi-set), meaning that each difference 
vector occurs with the same cardinality in either set; see |16| for an early class of examples 
in one dimension. 

A relatively simple homometric pair, realised as finite subsets Fi / F2 C 2 2 , was con- 
structed in [9] . One choice of the coordinates results in 



One can explicitly check that F\ — F\ = F2 — F^, including multiplicities. 

An appropriate generalisation to infinite point sets needs the concept of density. We call two 
infinite point sets homometric when their natural autocorrelation measures exist and coincide. 
Homometric point sets thus have the same density. Due to the volume averaging involved, 
two point sets related by adding or removing a point set of density are homometric. This 
is also true of point sets related by translation or inversion (but not, in general, by rotation). 

It is well known that two crystallographic (or fully periodic) point sets can only be homo- 
metric when they share the same lattice of periods. They are then mutually locally derivable 
(MLD) from each other [U[2], which also implies that the associated dynamical systems (un- 
der the translation action [15]) are topologically conjugate The corresponding question 
for mathematical quasicrystals (model sets without any periodicity) is more difficult, as we 
shall demonstrate by an example. 

2.1. Covariogram. There is an interesting connection between the homometry of model sets 
(with a Euclidean internal space) and the covariogram problem. For a non-empty, relatively 
compact subset K C which is assumed to be Riemann measurable, the function 



defined for all x £ M. d , is called the covariogram of K. The covariogram problem amounts 
to determine K from its covariogram cvg^-(x); compare [U [9]. This is sometimes also re- 
ferred to as Matheron's problem, which was originally formulated as the question whether 
the covariogram determines a convex body, among all convex bodies, up to translation and 
inversion; see [12 [131 [6] f° r details. Since cvg K (x) = cvg t+K (x) for any translation t G M. d and 
cvg^-(x) = cvg_^(x), the covariogram cvg^ can determine K at best up to translations and 
inversion. We call two non-empty, relatively compact, Riemann measurable sets K, K' C M d 
homometric when cvg^ = cvg^-,. 

Denoting the characteristic function of K by 1 K , the function cvg^ is given by the convo- 
lution 



(1) 



*i = { ffl > (0) .ffl.fi). (3) . (?) . (2) . (3) . (S . © . (3) . (') . S > (4) . (5) } . 
n ={(S). (?). (0). (i). fi). fi). (I), fi). (3), ©. Q. (3). ('). ©. (SI- 



(2) 




(3) 




Its Fourier transform 
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is an analytic, positive function that vanishes in the limit as \k\ — > oo. This relation is the 
reason why, if K is itself inversion symmetric in the sense that —K = t + K for a suitable 
translation t G R , the function 1 K (and hence K) can be reconstructed from the knowledge 
of cvg^- , up to translation and inversion [8] . 

If K is a convex poly tope in dimension d < 3, it is determined by cvg^; see [HI El EG] an d 
references therein. In general, however, the reconstruction of K from the knowledge of cvg^ 
is a difficult problem. An interesting example of two polyominoes with the same covariogram 
|3j follows from the point set pair of Eq. ([I]) by adding the unit square C = [— \ , ^] , so that 

(5) Px=Fx + C and P 2 = F 2 + C. 

Their covariograms are equal as a consequence of the homometry of the finite point sets Pi 
and F 2 , whence Pi and P 2 are homometric (as are also any translates of ±Pi and ±P2). The 
polyominoes Pi, P 2 and their joint covariogram are displayed and discussed in more detail in 
|3j. Let us mention in addition that the scaled polyominoes aPi and aP 2 are homometric to 
each other for any choice of a G R. 

2.2. Homometry of model sets. Let us now consider the situation of regular model sets 
A that are defined via a cut and project scheme [HI [2] with Euclidean internal space 

R d <-?— M. d x R m -^-> R m 
U U U dense 

(6) *(£) c > Wint{C ) 

L > V 

by 

(7) A = t+ X(W) = t + {x£ -k(C) j x* G W}. 

Here, £ is a lattice in R rf x R m , the window W C R m is a non-empty, relatively compact set 
with boundary of measure 0, and t G R d is an arbitrary translation. 

The autocorrelation j A of the corresponding Dirac comb 5a = Ylx&A $x exists and has the 
explicit form j A = YlzeA-A w ^ coefficients 

vol(Wn(W-z*)) 

(8) r,(z) = dens(A) V ^ {w) >± = dens(£) cvg w (z*), 

expressed in terms of the covariogram of the window W. Hence, two (Euclidean) model sets 
obtained from the same cut and project scheme are homometric if and only if the defining win- 
dows share the same covariogram. As a consequence, homometric model sets from the same 
cut and project scheme have the same diffraction measure. Conversely, kinematic diffraction 
cannot discriminate between homometric model sets. 

A planar example is obtained by using the homometric pair of polyominoes Pi and P 2 as 
windows for model sets in a cut and project scheme of type Q. For instance, we can use the 
Minkowski embedding C C 2 ~ R of L = Z[£s], where £s is a primitive 8th root of unity, 
and a *-map defined by a suitable algebraic conjugation. The two model sets A\ := A.(Pi) 
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and Ai := A(P2), with Pi and P 2 as defined in ([5]), are then homometric by construction 
(the relative position of the windows, which is irrelevant for homometry, maximises their 
intersection). The two model sets A\ and A<i are not locally indistinguishable, and differ in 
points of positive density. In particular, the difference sets A\ \A2 and A2 \A± are model sets 
themselves (but not homometric). 

The diffraction measure 7 is the same for both A\ and A2, and reads 

with intensity function I(y) = |ylj(y)| 2 derived from 

Ai(y) = dens(£ 8 ) l P .{-y). 

While the amplitudes depend on the window, their absolute squares do not. One can work out 
the explicit diffraction intensities; see [3] for details. The ratio of the (complex) amplitudes 
is given by 

A x (y) 1 + e 2"M>a + e *r»(Wi+2Wa) 1 - e^Vi 

( ' A 77a ~ 1 , o _-n-iO.li. 4-3l;-1 I \ ~~ 



A 2 (y) 1 + 2 e 7r *( 2 fi+ 3 s / 2) cos(7ry 2 ) e 27Tiy i + e 2 ™( y i +y z) + e~ 27Tiy z 

with y = (2/1,2/2)- This is a well-defined function on internal space M. 2 , with values in S , unless 
the denominator vanishes. The latter happens for y 2 G Z + {^, |} together with y 1 € Z. One 
can check that the ratio has no continuous extension to these points. Writing the ratio as 
exp(27rix(y)) j the phase function \ is n °t defined at these points. Moreover, as one can check 
explicitly, it does not satisfy the additivity property x{y + y') = x(y) J i~x(y') mod 1, wherefore 
the ratio fails to be a character on M 2 by violating both defining properties; compare [15J and 
references therein. Note that an analogous phenomenon already shows up in the comparison 
of Pi with —Pi, because these windows are not inversion symmetric up to translations. 

The choice of the windows Pi and P2 is special in the sense that A± and A2 turn out to be 
MLD, because the square C satisfies 

C = Pi n (-t + Pi) = P 2 n {-t + P 2 ) 

with the translation t = (4,5), and each window is now the union of 15 integral (and hence 
admissible) translates of C according to Eq. ([5]); see [1] for details. In this case, the associated 
dynamical systems |15j are again topologically conjugate. As mentioned above, the two 
windows may be scaled (by the same factor) without affecting their mutual homometry. For 
almost all choices of the scaling factor, one loses the MLD property of the corresponding 
model sets, because the finite reconstruction property [1] is lost. Nevertheless, the associated 
dynamical systems will always be metrically isomorphic (due to the Halmos-von Neumann 
theorem). It is an interesting open question whether they are still also topologically conjugate, 
which is a weaker equivalence notion than MLD. 

Independent of this conjugacy issues, our example illustrates that diffraction (hence auto- 
correlation) alone is generally insufficient to uniquely determine a regular model set. However, 
as discussed in [15], this ambiguity can be resolved with the knowledge of the 3-point corre- 
lations. This statement is immediate in our example (via the existence or non-existence of 
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certain patches), but holds in full generality for regular model sets; see [15] and references 
therein. 

3. Random Dirac combs 

The problem of reconstruction becomes even more involved in the case of mixed spectra. In 
this setting, a slight change in point of view is helpful to separate distinct spectral components. 
This is most easily achieved by considering weighted Dirac combs of point sets, with real 
(or even complex) weights. Below, generalising an example discussed in [10], we construct 
a family of one-dimensional homometric (weighted) point sets, based on the binary Rudin- 
Shapiro sequence, which cover the entire entropy range from to log (2), the maximal possible 
entropy for a binary system. This shows that, in general, it is not even possible to determine 
the degree of long-range order of the weighted point set from diffraction data. The same 
conclusion also holds for the diffraction of the associated unweighted point sets. 

3.1. Bernoulli versus Rudin-Shapiro. We start by re-considering the example of Ref. |10j . 
The first model is a Bernoulli system on Z, with the stochastic Dirac comb 

( 10 ) u B =^Y m 5 m , 

where iXm) m &z ls a f am ily of i.i.d. random variables that each take the values 1 and — 1, 
with probabilities p and 1 — p, where < p < 1. For the stochastic Dirac comb u> B , the 
autocorrelation measure 7 B and the diffraction measure 7^ almost surely exist and read 

7b = (2p-l) 2 S z + Ap(l-p)5 , 
7b= (2p-l) 2 5 z + 4p(l-p)A, 

where A denotes Lebesgue measure on R. Note that, in this stochastic situation, almost sure 
results are unavoidable. In particular, one has 7^ = A for p = ^ and 7^ = 5 Z for p = or 
p = 1. The choices p = and p = 1 correspond to the deterministic limiting cases oj b = ±5%, 
while p = \ describes a stochastic comb (coin tossing) with weights of average 0. 

The binary Rudin-Shapiro sequence is defined in two steps as follows |17j . We start from 
the substitution rule 

(12) a 1 — ^ ac, b 1— > dc, c 1— > ab, d 1— > da, 

on the four-letter alphabet A = {a, b, c, d}. We choose a bi-infinite fixed point (under the 
square of the above substitution, with seed ba) and apply the morphism tp : A — ► {±1} 
defined by ip(a) = ip(c) = 1 and tp(b) = tp(d) = —1, extended to A z . The autocorrelation and 
diffraction measures of the resulting binary Rudin-Shapiro chain S*rs are 

(13) 7 RS = 5 and 7^ = A. 

This is an example with a purely absolutely continuous diffraction, despite the fact that the 
Rudin-Shapiro chain is deterministic and has entropy 0. In particular, it agrees with the 
diffraction measure of the Bernoulli comb with p = ^, which has entropy log(2). 
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3.2. 'Bernoullisation'. It is possible to impose the influence of chance on the order of a 
deterministic system, and thus interpolate between deterministic and random systems. Here, 
we focus on binary sequences and modify them by an i.i.d. family of Bernoulli variables. 

Consider a bi-infinite binary sequence S € {±1} Z which we assume to be uniquely ergodic 
(in the sense that its hull under the action of the shift map is a uniquely ergodic dynamical 
system). Then, the corresponding Dirac comb uj s = Yliez^iSi possesses the (natural) au- 
tocorrelation 7 S = J2 m £Z r ls( m ) $m with autocorrelation coefficients rj s (m), where %(0) = 1 
by construction. 

Let (i^) igZ be an i.i.d. family of random variables that each take values +1 and — 1 with 
probabilities p and I — p. The 'Bernoullisation' of uj s is the random Dirac comb 

(14) ug^-^SiYiSi, 

i<=Z 

which emerges from u> s by independently changing the sign of each Si with probability 1— p. 
Setting Zi := S^X defines a new family of independent (though, in general, not identically 
distributed) random variables, with values in {±1}. Despite this modification, the auto- 
correlation 7 s . p of uj s . p almost surely exists and can be determined via its autocorrelation 
coefficients r] s . p (m) as follows. Since one always has rf s (0) = %(0) = 1, let m ^ and 
consider, for large N, the sum 

1 N 1 - 

i=~N (+,+) (-,-) (+,-) (-,+) 

which is split according to the value of (Si, Si- m ). Each of the four sums can be handled in 
the same way as for the Bernoulli comb, thus contributing (2p — l) 2 times the frequency of 
the corresponding sign pair. Observing that the overall signs are the products SiSi- m , it is 
clear that, as N — > oo, one obtains (almost surely) 

(16) 7 ls . p (m) = (2p-l) 2 Vs (m) 

for all Thus, the autocorrelation ^ s . p of uJ s . p almost surely exists and is given by 

(17) 1S;p = (2p-l) 2 ls + Ap(l-p)5 , 

where 7 5 is the unique autocorrelation of u> s . 

Consider now the Bernoullisation (with parameter p) of the binary Rudin-Shapiro sequence 
with random Dirac comb u; RS . Its autocorrelation measure almost surely exists and reads 
7 RS . p = <5 , independently of p. This means that the random Dirac combs w RS . p , even for 
different values of p, are almost surely homometric, and share the purely absolutely continuous 
diffraction measure 7^^, = A. 

3.3. Extension to higher dimension. Our above discussion has an immediate extension 
to Euclidean space of arbitrary dimension d. Consider d complex- valued sequences (U^) ieZ 
and define the weighted Dirac comb 

^ = [/(\ 8 ...^(\ = ^(nc/(;))t 
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Here, x = . . . ,x d ), and the weights (on Z, d ) are products of d elements of the individual 
sequences. Assuming that the natural autocorrelations of the individual sequences exist, the 
relevant observation is that the resulting autocorrelation of coy (and hence also the corre- 
sponding diffraction measure) factorises accordingly. 

Each U® may be chosen as a member of our previous one-parameter family of (|14p . in 
particular as C/w = Srs f° r an £■ This results in a deterministic weighted Dirac comb on 
Z d with diffraction measure 7 = A, where A now denotes Lebesgue measure on M. d . This 
represents a system of entropy 0. 

The analogue of the Bernoullisation of (fl~4"j) in d-space, with an i.i.d. family (%:) x€ ^d, then 
results in an isospectral one-parameter family of random Dirac combs on which can realise 
any entropy between and log(2). Similarly, using 7L d ~ Z fc x TL d ~ k for some < k < d and 
restricting the Bernoullisation to Z fc , one obtains isospectral families with arbitrary entropy 
of rank k between and log(2). This indicates that the variety of homometric structures 
grows with dimension. 

4. Conclusions 

The homometry problem for regular model sets in dimensions d < 3 appears to have a 
unique solution if one may assume that the window is convex [HE]. However, this favourable 
situation is not always met in real quasicrystals. Our example, with non-convex windows, 
illustrates the existence of distinct homometric structures. Since homometric model sets are 
always metrically isomorphic, but not necessarily MLD, it remains an interesting question 
how they are related as topological dynamical systems. 

For the case of spectra with a continuous component, the Bernoullisation approach can 
explore the full entropy range: the Bernoulli case (with p = i) has entropy log 2, the maximal 
value for a binary system, while Rudin-Shapiro has entropy 0, and the parameter p interpolates 
continuously between the two limiting cases. This gives an indication of how degenerate the 
inverse problem really is. Unless additional information is available, one possible strategy 
to proceed would employ an optimisation approach, for instance by choosing the structure 
which maximises the entropy, which singles out the Bernoulli comb here. 
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